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ABSTRACT 

We compute the luminosity function (LF) and the formation rate of fong gamma ray 
bursts (GRBs) by fitting the observed differential peak flux distribution obtained by 
the BATSE experiment in two different scenarios: i) the GRB luminosity evolves with 
redshift and ii) GRBs form preferentially in low-metallicity environments. In both 
cases, model predictions are consistent with the Swift number counts and with the 
number of detections at z > 2.5 and z > 3.5. To discriminate between the two evolu- 
tionary scenarios, we compare the model results with the number of luminous bursts 
(i.e. with isotropic peak luminosity in excess of 10^"^ erg s~^) detected by Swift in 
its first three years of mission. Our sample conservatively contains only bursts with 
good redshift determination and measured peak energy. We find that pure luminosity 
evolution models can account for the number of sure identifications. In the case of a 
pure density evolution scenario, models with Zth > 0.3 Zq are ruled out with high 
confidence. For lower metallicity thresholds, the model results are still statistically 
consistent with available lower limits. However, many factors can increase the discrep- 
ancy between model results and data, indicating that some luminosity evolution in 
the GRB LF may be needed also for such low values of Zth. Finally, using these new 
constraints, we derive robust upper limits on the bright-end of the GRB LF, showing 
that this cannot be steeper than 2.6. 
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1 INTRODUCTION 

Gamma Ray Bursts (GRBs) are powerful flashes of high- 
, energy photons occurring at an average rate of a few per day 
throughout the universe. Their luminosity is such that they 
can be detected up to very high redshift (the current record 
is GRB 080913 at z = 6.7; Greiner et al. 2008). The energy 
source of a GRB is believed to be associated to the collapse 
of the core of a massive star in the case of long-duration 
GRBs, and due to merger- or accretion-induced collapse 
for the short-hard class of GRBs (see Meszaros 2006 for a 
recent review). In this paper, we limit our analysis to the 
class of long-duration GRBs. 

The knowledge of GRBs has enormously benefited from 
the observations of the Swift satellite (Gehrels et al. 2004) . 
Although the current sample of GRBs with known redshift 
is still too poor to allow a direct measure of the GRB lu- 
minosity function (LF), important constraints on the cos- 
mic evolution of these sources can be set on the basis of 
recent Swift results. In particular, Salvaterra & Chincarini 
(2007, herethereafter SC07) showed that models in which 
GRBs are unbiased tracer of cosmic star formation and are 



characterised by a constant LF are robustly ruled out by 
the number of GRB detections aX z > 2.5 and z > 3.5. 
Similar conclusions were reached recently by other studies 
(e.g. Guetta, Piran & Waxman 2005; Daigne et al. 2006; 
Gen & Fang 2007; Kistler et al. 2008). Moreover, they have 
shown that Swift data can be reproduced assuming lumi- 
nosity evolution of the GRB LF (see also Lloyd-Ronning, 
Fryer & Ramires-Ruiz 2002; Wei & Gao 2003; Daigne et al. 
2006) and/or that GRBs form preferentially in low metallic- 
ity environments (see also Natarajan et al. 2005; Langer & 
Norman 2006; Li 2007; Gen & Fang 2007; Lapi et al. 2008). 

In this Letter, we derive the formation efficiency and 
the free parameters describing the GRB LF by fltting the 
differential peak flux distribution of BATSE GRBs in these 
two scenarios. We then obtain new and tighter constraints 
on the cosmic evolution of long GRBs by comparing different 
models against the number of luminous (i.e. with isotropic 
peak luminosities L > 10^"^ erg s~^) GRBs detected by Swift. 
We also consider models with joint luminosity and density 
evolution providing a robust upper limit on the steepness of 
the bright-end of the GRB LF. 

This Letter is organised as follows. In Section 2, we 
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briefly describe the different models and the main equations 
used in the calculation of the LF, and in Section 3, we com- 
pare model results against Swift data. Finally, we summarise 
our findings in Section 4. 



2 MODEL DESCRIPTION 

The observed photon flux, P, in the energy band 

E < -Bmax, emitted by an isotropically radiating source at 

redshift z is 



P = 



(1) 



where S{E) is the differential rest-frame photon luminos- 
ity of the source, and dL(z) is the luminosity distance. 
To describe the typical burst spectrum we adopt the func- 
tional form proposed by Band et al. (1993), i.e. a smoothly 
broken power-law with a low-energy spectral index a, a 
high-energy spectral index (5, and a break energy Eh = 
(a - 0)Ep/(2 + a), with a = -1 and 3 = -2.25 (Prcccc ct 
al. 2000, Kaneko et al. 2006). The spectrum normalisation 
is obtained by imposing that the isotropic-equivalent peak 
luminosity is L = jl'^^^^'^^ ES(E)dE. In order to broadly 
estimate the peak energy of the spectrum, Ep, for a given 
L, we assumed the validity of the correlation between Ep 
and L (Yonetoku et al. 2004, Ghirlanda et al. 2005), which 



is basically a diffenent expression of the Ep 
(Amati et al. 2002, Amati 2006). 



: 337 keV 



2 X lO-''^ erg s" 
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Although the above correlation has an appreciable scatter, 
we will show that this does not affect our results. 

Given a normalised GRB LF, 4>{L), the observed rate 
of bursts with peak flux between Pi and P2 is 
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where dV{z)/dz = 47rcd|,(2)/[i/(a)(l -|- z)^] is the comov- 
ing volume element*, and H{z) = Ho[ilM{l + z)^ + CIa + 
(1 — ~ 0,a){1 + z'f]^^"'^. ASls is the solid angle covered 
on the sky by the survey, and the factor {1 + z)~^ accounts 
for cosmological time dilation. Finally, ^'grb(2) is the co- 
moving burst formation rate. In this work, we model the 
GRB LF with a power law with an exponential cut-off at 
low luminosities: 



4>{L) cx ( -— 
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-^cut 

IT 



(4) 



SC07 have shown that models in which GRB form pro- 
portionately to the star formation rate (SFR) and are de- 
scribed by a LF constant in redshift are robustly ruled out 
by the number of GRB with sure detection at 2; > 2.5 and 

* We adopted the 'concordance' model values for the cosmologi- 
cal parameters: h = 0.7, Clm = 0.3, and Qa = 0.7. 



Pure Luminosity Evolution Models 
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1.5 


1.07±0.11 


0.66±0.21 


2.16±0.09 


0.95 


2.0 


l.OlitO.lO 


0.36±0.08 


2.08±0.06 


0.83 


2.5 


0.94±0.09 


0.20±0.04 


2.03±0.05 


0.92 


3.0 


0.94±0.09 


0.10±0.02 


1.99±0.04 


0.84 


Pure Density Evolution Models 




feGRB/(10-8M-l) 


Lo/(105l erg s-i) 




Xr 


0.1 


11.38±1.60 


10.34±4.10 


2.51±0.22 


0.81 


0.2 


4.45±0.68 


8.40±3.58 


2.48±0.22 


0.78 


0.3 


2.80±0.45 


7.40±3.11 


2.50±0.22 


0.80 



Table 1. Best— fit parameters for difi'erent models: top panel for 
pure luminosity evolution and bottom panel for pure density evo- 
lution. Errors are at Itr level. 



z > 3.5 during the first two years of Swift mission. Thus, 
GRBs should have oxporiencod some kind of evolution, be- 
ing more luminous or more common in the past. Therefore 
wo consider here two families of models: (i) luminosity evo- 
lution models, where the cut-off luminosity in the GRB LF 
varies as Lcut = Lo(l-|-2;)'' and (ii) density evolution models, 
where GRBs form preferentially in galaxies with metallicity 
below a given threshold Zth- In the first case, the GRB for- 
mation rate is simply proportional to the global SFR, i.e. 
^^GRBiz) = kGRB'^*{z). We use here the recent determi- 
nation of the SFR obtained by Hopkins & Beacom (2006), 
slightly modified to match the observed decline of the SFR 
with {l + z)~^'^ at 2; > 5 suggested by recent deep-field data 
(Stark et al. 2006). For the density evolution models, the 
GRB formation rate is obtained by convolving the observed 
SFR with the fraction of galaxies at redshift z with metal- 
licity below Zth using the expression computed by Langcr 
&: Norman (2006). In this scenario, the cut-off luminosity is 
assumed to be constant in redshift, i.e. Lcut = const = Lq. 

Furthermore, we consider a third family of models in 
which both effects are presented: GRBs form preferentially 
in environments with Z < Zth and are characterised by an 
evolving LF. 



3 LUMINOSITY vs DENSITY EVOLUTION 

The free parameters in our model arc the GRB formation 
efficiency fccRB, the cut-off luminosity at 2; = 0, Lo, and 
the LF power index ^. We optimised the value of these pa- 
rameters by minimisation over the observed differential 
number counts in the 50-300 keV band of BATSE. We use 
here the results by Stern et al. (2000) , who considered both 
triggered and non-triggered bursts and also corrected the 
distribution taking into account the BATSE detector effi- 
ciency. The best-fit parameters are reported in Table 1. As 
already pointed out by SC07, it is always possible to find a 
good agreement between models and data. Moreover, it is 
possible to reproduce also the differential peak fiux count 
distribution in the 15-150 keV Swift band using the same 
GRB LF and formation efficiency obtained by fitting the 
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Figure 1. Cumulative number of GRBs a,t z > 2.5 (left panels) 
and z > 3.5 (right panels) as a function of the photon flux P. 
The no-evolution case is shown with the solid thin line. Top panels 
refer to luminosity evolution models: solid line is for 5 = 3, dotted 
line for 5 = 2.5, short-dashed line for 5 = 2, and long-dashed for 
5 = 1.5. Bottom panels refer to density evolution models: solid 
line is for Zth = 0.1 Zq, dotted line is for Z^f^ = 0.2 Zq, and 
short-dashed line is for Z^^ = 0.3 Zq. The number of bursts 
detected by Swift in three years is shown as solid histogram. Note 
that the observed detections arc lower limits, since many high-2 
GRBs can be missed by optical follow-up searches. A field of view 
of 1.4 sr for Swift/B AT is adopted. 

BATSE data. Consistently with SC07, we find that the num- 
ber of GRBs confirmed at z > 2.5 and z > 3.5 in the three 
years of Swift mission requires some kind of evolution. The 
results arc shown in Fig. 1 for different evolution models. As 
a comparison, we show also the result for the no-cvolution 
model with the solid thin line. The chance probability as- 
sociated to the no-evolution model is found to be less than 
10"'' ensuring that this kind of model can bo discarded at a 
very high confidence level and indicating the need of some 
kind of evolution to explain Smft high-2; detections (see also 
SC07). 

In this work wc liiglilight a new toll to discriminate be- 
tween these two evolution scenarios by computing the num- 
ber of luminous GRBs, i.e. bursts with isotropic peak lumi- 
nosity L > 10^* erg s"^ in the 1-10000 keV band. The model 
predictions are obtained by 

where Lraax{z) — ma,x{L, Loa{z)} and Loa{z) is the mini- 
mum luminosity of a burst exploding at rcdsliift z able to 
trigger Swift, i.e. P{Lo.4,z) = 0.4 ph s~^ cm~^ in the 15-150 
keV BAT bandt. 

We compare model predictions with the number of 

^ We have assumed here a trigger threshold of 0.4 ph cm~-^ 
for which the Swift/BAT sample is complete: indeed, we find that 
the observed differential peak flux distribution of Swift GRBs 



L [erg s"'] L [erg s"'] 

Figure 2. The cumulative number of luminous GRBs detected by 
Swift in three years of operations is shown with the histogram as 
a function of the isotropic equivalent peak luminosity, L. Shaded 
area takes into account the errors on the determination of L. 
Note that the data are to be considered as lower limits of the 
real number of Swift detections, since only bursts with measured 
redshift and peak energy are included in our sample. Pure lu- 
minosity (density) evolution models are shown in the left (right) 
panel. Lines as in Fig. 1. A field of view of 1.4 sr for Swift/B AT 
is adopted. 

bright bursts detected by Swift in three years of mission. 
Conservatively, our data sample contains only bursts with 
a good redshift measurement and whose peak energy was 
measured or well constrained by Swift itself or other satel- 
Utes (such as HETE-2 or Konus-Wind). We find nine GRBs^ 
detected by Swift with L > 10'"''^ erg s~^ in three years of 
mission. We want to stress here that this number represents 
a conservative lower limit on the real number of bright GRBs 
detected, since some luminous bursts without z and/or Ep 
can be present in the Swift catalogue. In particular, we note 
that a good redshift determination is obtained for ~ 1/3 of 
Swtft bursts and for only a fraction of these we have a well 
constrained Ep. The cumulative distribution of the known 
bright bursts is shown in Fig. 2. The shaded area takes into 
account errors on the determination of L. The model results 
for the pure luminosity (density) evolution models are plot- 



below this threshold is less populated than what expected from 
fitting the BATSE data (see also SC07). 

■t Six out of the nine GRBs included in our sample are reported 
in Rossi et al. (2008): GRB 050401, GRB 050603, GRB 060927, 
GRB 061007, GRB 061121, and GRB 071020. In this work we 

use the peak luminosity computed on a 1 s timcscalc. Wc add to 
these other three bursts that arc not present in Rossi ct al. (2008) 
compilation: GRB 050505 with L = 1.6 ± 1.0 X 10"'=^ erg s^l, 
GRB 060210 with L = 4.6 ±2.8 x lO^^ erg s"!, and GRB 060124 
with L = 1.1 ± 0.1 X 10^53 erg s"^ (Romano et al. 2006). We 
note that the peak fluxes of all of the nine bursts considered 
in our sample are well above the assumed trigger threshold of 
Swift/BAT. 
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ted in the left (right) panel of Fig. 2. We have also checked 
that our findings do not depend on the assumed L — Ep cor- 
relation: considering a mean break energy (as done in SC07) 
for all bursts does not change significantly our results. 

For the pure luminosity evolution scenario, all models 
here considered predict a large number of bright bursts to 
be detected by Swift. Indeed, they can easily account for the 
observed number of bursts with L > 10^"' erg s~^. On the 
other hand, models in which GRB formation is confined in 
low-met allicity environments seems to fall short to account 
for the observed bright GRBs for Zth > 0.1. In particular, for 
Zth = 0.3 Zq , as required by collapsar models (MacFadyen 
& Woosley 1999; Izzard, Ramirez- Ruiz & Tout 2004), only 
~ 6 bursts with L > 10'''^ erg s^'^ should have been detected 
in three years of observations. 

In order to test the statistical significance of our find- 
ings, wc explore the parameter space around the best fit 
parameters looking for triple (/cgrb, -^o) compatible with 
the sure bright burst identifications. Among these, we choose 
the triple that give the best agreement with the differential 
peak flux distribution of BATSE GRBs. Then, the null hy- 
pothesis test gives us the confidence level at which we can 
discard the considered model. For Zth ~ 0.3 Z©, wc find 
that null hypothesis is satisfied being, in the best case, the 
chance probability of ~ 0.22. However, since we are deal- 
ing with strong lower limit on the real number of bright 
burst detections, we have also to consider the case in which 
some luminous bursts are hidden among the GRBs missing 
redshift and/or Ep measurement. We find that the null hy- 
pothesis probability decreases rapidly with the number of 
bright burst detections and already for 10 bright bursts it 
drops down to the percent level (see Fig. 3). So, although 
the Zth = 0.3 Zq scenario cannot formally be discarded by 
available data, some degree of luminosity evolution in the 
GRB LF is suggested if just a few more bright bursts had 
to be added to our measured sample. For Zth =0.1 (0.2) 
Zq, the models can account up to 16 (12) bright bursts in 
the whole Swift data sample, i.e. just 7 (3) 'hidden' l)riglit 
bursts (Fig. 3). For the pure evolution scenario, we find that 
even a relatively large population of bright bursts hidden in 
the Swift catalogue of bursts without good redshift and/or 
Ep measurement, can still be accounted for. 

To complete our analysis, we then consider models in 
which both luminosity and density evolution are present. 
For any value of the threshold metallicity Zth, we derive the 
minimum luminosity evolution required to match the num- 
ber of known GRBs with L > 10^^ erg in the 3-year 
catalogue (top panel of Fig. 4). The error bars correspond 
to a null hypothesis probability of 1% that the model re- 
produces the differential peak flux distribution of BATSE 
GRBs. We find that for Zth > 0.3 Z©, the available con- 
straints require luminosity evolution in the GRB LF. Below 
this threshold, models charaterised by a constant LF can ac- 
count for the number of known luminous bursts considering 
the error bars. As already pointed out, the existence of a few 
bright bursts hidden among GBRs without a sure redshift 
and /or Ep measure (more than 2/3 of the whole Swift cat- 
alogue) , would imply the need of luminosity evolution even 
for lower value of Zth- Moreover, note here that the exis- 
tence of a distinct population of long low-luminosity GRBs 
at low-z (see e.g. Guetta & Delia Valle 2007) strengthens our 
conclusions. In the form of the adopted LF (see eq. (4)), we 
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Figure 3. Chance probability to reproduce the differential peak 
flux distribution of BATSE GRBs given a model that predicts 
N bursts with L > 10^3 erg s ^ in 3 years of Swift operations. 
Different pure density evolution models are plotted with different 
lines: dashed line for Z^h = 0.3 Zq, dotted line for Z^h = 0.2 Zq, 
and solid line for Zth = 0.1 Zq. The probability rapidly decreases 
with the number of bright bursts so that, even a few more bright 
bursts to be added to our conservative sample would require some 
degree of luminosity evolution. 



do not include this population in our calculations. Should 
this population be statistically significant and present at 
all redshifts, the faint end of the LF would be more pop- 
ulated and all model predictions would be shifted towards 
lower values, increasing the discrepancy between model re- 
sults and data. Since for pure density CA'olution models the 
LF cut-off luminosity is considerably larger than for pure 
luminosity evolution models, the former models would be 
more severely affected by the existence of a largo population 
of undcrluminous bursts at any redshift. In conclusion, al- 
though wc cannot rule out at a high confidence level pure 
density evolution models with Zth < 0.3 Z©, available data 
suggest some luminosity evolution in the GRB LF with red- 
shift. 

Finally, we constrain the steepness of the bright-end of 
the GRB LF (bottom panel of Fig. 4). Very steep LF are 
robustly ruled out in every model considered. We find that 
the maximum value of the index ^ have is < 2.6. For the pure 
luminosity evolution models, wc find ^ < 2.2. As wc already 
pointed out, these limits could be further constrained to 
explain current Swift data due to the possible existence of 
bright bursts missed by our conservative selection and of 
a relatively large population of faint bursts not considered 
here. 



4 CONCLUSION 

We have computed the luminosity function and the forma- 
tion rate of long GRBs by fitting the observed differential 
peak flux distribution obtained by the BATSE satellite in 
two different scenarios: i) the GRB luminosity evolves with 
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0.6 0.8 

Figure 4. Top panel: lower limit on the luminosity evolution 
parameter, 5, for different threshold metallicity Zth- Error bars 
correspond to a null hypothesis probability, that the model re- 
produces the differential peak flux distribution of BATSE GRBs, 
equal to 1%. Bottom panel: corresponding upper limit on the 
bright-end power index of the GRB LF, ^. 



redshift and ii) GRBs form preferentially in low-metallicity 
environments. In both cases, model predictions are consis- 
tent with the Swift number counts and with the number of 
detections at z > 2.5 and z > 3.5. To discriminate between 
the two evolutionary scenarios, wc compared the model re- 
sults against the number of luminous bursts, i.e. with peak 
luminosity in excess of 10®^ erg s~^, detected by Swift. Con- 
servatively, our data sample contains only bursts with good 
redshift determination and measured peak energy. We find 
that models in which GRBs are characterised by a constant 
LF (i.e. for pure density evolution models) are disfavoured 
as they underpredict the number of luminous GRBs. Swift 
data can be explained assuming that the GRB luminosity 
evolves with redshift. Although we cannot discard pure den- 
sity evolution models with Zth < 0.3 on the basis of the 
current sample, the existence of a few bright GRBs missed 
by our conservative selection criteria and/or of a relatively 
large population of faint GRBs would require some luminos- 
ity evolution in the GRB LF even for such low values of Zf^. 
On the other hand, pure luminosity evolution scenarios can 
account more easily for a largo number of burst detections 
with L > 10''"' erg s^^. In this work, we derive lower limits 
to the luminosity evolution of the GRB LF with redshift for 
different values of the metallicity threshold for the GRB for- 
mation. Moreover, we use these constraints to set a robust 
upper limit on the bright -end of GRB LF. Wo find that the 
number of bright GRBs detected by Swift implies that this 
cannot be very steep: ^ < 2.2 (pure luminosity evolution) 
and i < 2.6 (for Ztn < 0.3 Zq). 

In conclusion, we find that available Swift observations 
point towards a scenario whore GRBs wore more luminous in 
the past. Although current data sample of bright GRBs with 
good redshift and Ep determination is still very poor, our 
findings show that these data can be used to set important 



constraints on the cosmic evolution of GRBs and on the 
steepness of their LF. 

Finally, the new constraints on the GRB LF allow us to 
derive robust lower limits on the number of bursts detectable 
by Swift at very high redshift. Assuming a trigger threshold 
Piim = 0.4 ph s~^ cm~^, at least ~ 5 — 10% of all detected 
GRBs should lie at ^: > 5, where the lower (upper) value 
refers to a pure luminosity evolution (pure density evolution 
with Zth = 0.1 Zq) model. Among these, > 1 — 3 GRB yr~^ 
should be detected at 2: > 6. These lower bounds double 
by lowering the Swift trigger threshold by a factor of two 
(Salvaterra et al. 2008). These results are consistent with 
the lower limits on the number of high-2: detections obtained 
by Salvaterra et al. (2007, 2008) using redshift distribution 
constraints. 
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